Dissociation of lp quarkonium states in a hot QCD medium 
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We extend the analysis of a very recent work (Phys. Rev. C 80, 025210 (2009)) to study the 
dissociation phenomenon of lp states of the charmonium and bottomonium spectra (xc and Xb) in a 
hot QCD medium. This study employed a medium modified heavy quark potential which is obtained 
by incorporating both perturbative and non-perturbative medium effects encoded in the dielectric 
function to the full Cornell potential. The medium modified potential has a quite different form (a 
long range Coulomb tail in addition to the usual Yukawa term) compared to the usual picture of 
Debye screening. We further study the flavor dependence of their binding energies and dissociation 
temperatures by employing the perturbative, non-perturbative, and the lattice parametrized form 
of the Debye masses. These results are consistent with the predictions of the current theoretical 
works. 
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I. INTRODUCTION 

One of the amazing discoveries of experimental mea- 
surements at RHIC is the surprising amount of both 
radial [l[ and elliptic flow 0] exhibited by the outgo- 
ing hadrons. Theoretical calculations cannot generate 
sufficient flow to explain the observations unless par- 
tonic cross sections are artificially enhanced by more 
than an order of magnitude over perturbative QCD 
predictions Q . Thus the matter created in these collisions 
is strongly interacting, unlike the type of weakly interact- 
ing quark-gluon plasma expected to occur at very high 
temperatures on the basis of asymptotic freedom [4j . The 
behavior of the heavy quarkonium states in hot strongly 
interacting matter was proposed as test of its confinement 
status, since a sufficiently hot dcconfincd medium will 
dissolve any binding between the quark-antiquark pair 
[j| . Another possibility of dissociation of certain quarko- 
nium states (sub-threshold states at T = 0) is the decay 
into open charm (beauty) mesons due to in-medium mod- 
ification of quarkonia and heavy- light meson masses Q. 

Many attempts have been made to understand the dis- 
sociation phenomenon of QQ states in the dcconfincd 
medium, using either lattice calculations of quarkonium 
spectral functions @, H, & EE EH or non-relativistic cal- 
culations based upon some effective potential [HI, EE EE 
EE EE E3- These two approaches show poor match- 
ing between their predictions because of the uncertain- 
ties coming from a variety of sources. None of the ap- 
proaches give a complete framework to study the proper- 
ties of quarkonia states at finite temperature. However, 
some degree of qualitative agreement had been achieved 
for the S-wave correlators. The finding was somehow am- 
biguous for the P-wave correlators and the temperature 
dependence of the potential model was even qualitatively 



different from the lattice one. Refinement in the com- 
putations of the spectral functions have recently been 
done by incorporating the zero modes both in the S- and 
P-channels [IE, E^- It was shown that, these contribu- 
tions cure most of the previously observed discrepancies 
with lattice calculations. This supports the use of poten- 
tial models at finite temperature as an important tool to 
complement lattice studies. 

The production of J/ip and T mesons in hadronic re- 
actions occurs in part through production of higher ex- 
cited cc (or bb) states and their decay into quarkonia 
ground state. Since the lifetime of different sub-threshold 
quarkonium states is much larger than the typical life- 
time of the medium which may be produced in nucleus- 
nucleus collisions; their decay occurs almost completely 
outside the produced medium. This means that the pro- 
duced medium can be probed not only by the ground 
state quarkonium but also by different excited quarko- 
nium states. Since, different quarkonium states have dif- 
ferent sizes (binding energies), one expects that higher 
excited states will dissolve at smaller temperature as 
compared to the smaller and more tightly bound ground 
states. These facts may lead to a sequential suppression 
pattern in J/ip and T yield in nucleus-nucleus collision as 
the function of the energy density. So, if one wants to in- 
terpret the J/ip suppression pattern observed in nuclear 
collisions at CERN SPS and RHIC, as a signature of the 
formation of the QGP, one requires a right understand- 
ing of the dissociation of Xc and Xb in the QGP medium. 
This is due to the fact that a significant fraction (~ 30%) 
of the 3 1 tjj yield observed in the collisions is produced by 
Xc decays [2fJ, |2l|, |22j . The J/ ip yield could show a signif- 
icant suppression even if the energy density of the system 
is not enough to melt directly produced J/ ?/> but it is suf- 
ficient to melt the higher resonance states because they 
are loosely bound compared to the ground state J/ifj. 
This motivates the special attention to the excited states 
Xc and Xb- 

In the studies of the bulk properties of the QCD plasma 
phase [IE HE HE HH , deviations from perturbative cal- 
culations were found at temperatures much larger than 
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the deconfinement temperature. This calls for quanti- 
tative non-perturbative calculations. The phase transi- 
tion in full QCD appears as a crossover rather than a 
'true' phase transition with related singularities in ther- 
modynamic observables (in the high-temperature and 
low density regime) [2lJ ■ Therefore, it is not reasonable 
to assume that the string-tension vanishes abruptly at 
or above T c and one should study its effect on the be- 
havior of quarkonia even above the deconfinement tem- 
perature. This issue, usually overlooked in the litera- 
ture, was certainly worth to investigate. This is exactly 
what we have done in our recent work [28l . [2^ | where we 
have obtained the medium-modified form of the heavy 
quark potential by correcting the full Cornell potential 
(linear plus Coulomb), not only its Coulomb part alone 
as usually done in the literature, with a dielectric func- 
tion encoding the effects of the deconfined medium. We 
found that this approach led to a long-range Coulomb 
potential with an (reduced) effective charge [28| in addi- 
tion to the usual Debye-screened form employed in most 
of the literature. With this effective potential, we inves- 
tigated the effects of perturbative and non-perturbative 
contributions to the Debye mass on the dissociation of 
quarkonium states. We subsequently used this study to 
determine the binding energies and the dissociation tem- 
peratures of the ground and the first excited states of 
charmonium and bottomonium spectra. 

However, our starting potential (Cornell) at T = 
has no terms to account the spin-dependence forces in 
QCD (jftj, so the medium-modified potential [28| also has 
no spin-dependent terms. As a consequence, Schrodingcr 
equation with the above medium-modified potential gives 
the same energy eigenvalues for the first excited states 
ip 1 and \c making them degenerate. This is certainly not 
desirable since their masses are not the same (in fact, 
mass of ifj' is slightly higher than \c whereas the latter 
state is more tightly bound than the former). Therefore 
the determination of the binding energies of lp states, 
viz., Xc, Xb and their dissociation temperatures like the 
ground and first excited states is not directly possible 
as had been done in our earlier work by employing the 
medium modified potential [28| . The principal quantum 
number (n) of ip' and Xc are same but their spin quantum 
number and as well as their total angular momentum are 
not the same. So, their quantum states should be denoted 
by all four quantum numbers (nlsj) and the difference in 
their binding energies (or in their total masses) should 
be originated from a spin-dependent correction terms. 

We have done this job in a two fold way. First, we 
have determined the binding en ergy for ip' by employing 
the medium-modified potential [28[ into the Schrodingcr 
equation. Then we obtain the binding energy for Xc 
by adding the correction terms to the binding energy 
of ip' . In our analysis, correction terms will be obtained 
by adopting a variational treatment of the relativistic 
two-fermion bound-states in quantum electrodynamics 
(QED) [33, taking into account the spin-dependent 
terms for the corresponding quantum numbers of ip' and 



Xc states. In this endeavor, coupled integral equations 
for a relativistic two-fermion system are derived varia- 
tionally within the Hamiltonian formalism of QED using 
an improved ansatz that is sensitive to all terms in the 
Hamiltonian 32 1 . 

The paper is organized as follows. In Sec. II, we re- 
view the work on the medium modified Cornell potential 
and dissociation of Is and 2s states of charmonium and 
bottomonium spectra. In Sec. Ill, we discuss how to de- 
termine the binding energies of Xc and Xb- in Sec. IV, 
we study the melting of Xc and Xb in the QGP medium 
and determine their dissociation temperatures. Finally, 
we conclude in Sec.V. 



II. IN-MEDIUM MODIFICATIONS TO 
HEAVY-QUARK POTENTIAL 

The interaction potential between a heavy quark and 
antiquark gets modified in the presence of a medium and 
it plays a vital role in understanding the fate of quark- 
antiquark bound states in the QGP medium. This issue 
has well been studied and several excellent reviews ex- 
ist [HI, Hl[ which dwell both on the phenomenology as 
well as on the lattice QCD. In these studies, they as- 
sumed the melting of the string motivated by the fact 
that there is a phase transition from a hadronic matter 
to a QGP phase. As a consequence they modified the 
Coulomb part of the potential only so they used a much 
simpler form (screened Coulomb) of the medium modi- 
fied potential in the deconfined phase. But recent lattice 
results indicates that there is no genuine phase transi- 
tion at vanishing baryon density, it is rather a cross-over, 
so there is no reason to assume the melting of string at 
the deconfinement temperature. We have addressed this 
issue in our recent work [28[ where we developed an effec- 
tive potential once one corrects the full Cornell potential 
with a dielectric function embodying medium effects. We 
recall the basic details which are relevant for the present 
demonstration. 

Usually, in finite-temperature QFT, medium modifi- 
cation enters in the Fourier transform of heavy quark 
potential as 



V(k) 



V(k) 
e(k) 



(D 



where e(k) is the dielectric permittivity given in terms 
of the static limit of the longitudinal part of gluon self- 
energy dH]: 



e(k) = 1 



n L (o,fc,T)\ _ 



k' 2 



= 1 



■#)■ (2) 



The quantity V(k) in Eq.|T|) is the Fourier transform 
(FT) of the Cornell potential. The evaluation of the FT 
of the Cornell potential is not so straightforward and can 
be done by assuming r- as distribution (r — > r cxp(— jr)). 
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After the evaluation of FT we let 7 tends to zero. Now 
the FT of the full Cornell potential can be written as 



V(fc) 



4a 



(3) 



Substituting Eqs.([2]) and (J3]) into ([I]) and then evaluat- 
ing its inverse FT one obtains the r-dependence of the 
medium modified potential [29| : 



V(r,T) 



2a \ exp (—mor) 
— ~ a 



2a 



' D 



2a 

m D 
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(4) 



This potential has a long range Coulombic tail in addi- 
tion to the standard Yukawa term. The constant terms 
are introduced to yield the correct limit of V(r,T) as 
T — ► (it should reduce to the Cornell form). Such 
terms could arise naturally from the basic computations 
of real time static potential in hot QCD[37| and from the 
real and imaginary time correlators in a thermal QCD 
medium [3^]. 

It is worth to note that the potential in a hot QCD 
medium is not the same as the lattice parametrized 
heavy quark free-energy in the deconfined phase which 
is basically a screened Coulomb [H, EFj] because one- 
dimensional Fourier transform of the Cornell potential 
in the medium yields the similar form as used in the lat- 
tice QCD to study the quarkonium properties which as- 
sumes the one-dimensional color flux tube structure [4l[ . 
However, at finite temperature that may not be the 
case since the flux tube structure may expand in more 
dimensions 39]. Therefore, it is better to consider the 
three-dimensional form of the medium modified Cornell 
potential which have been done exactly in the present 
work. We have compared our in-medium potential with 
the color-singlet free-energy [42j extracted from the lat- 
tice data and found that it agrees with the lattice re- 
sults except from the non-perturbative result of the De- 
bye masses. 

However, if we neglect the finite range terms viz. 
Yukawa term in the limit r > > 1 /tud and for large values 
of temperatures the product am o will be much greater 
than 2a /mo, then it leads to an analytically solvable 
Coulomb potential: 



V(r,T) 



2a 



amjj 



(5) 



We employed this medium-modified effective potential 
to study the binding energies and the dissociation tem- 
peratures for the ground and first excited states of cc 
and bb spectroscopy. However, to see the effects of the 
finite-range terms, we solve the Schrodingcr equation nu- 
merically with the full effective potential(|l]) and found 
that the dissociation temperatures was changed by ~ 
10% (2q . Therefore the approximated form ([5|) have 
a dominant role in deciding the fate of these states in 



the hot QCD medium. Let us now proceed to the de- 
termination of the binding energies and the dissociation 
temperatures for \c and Xb states in Sec(s).III and IV, 
respectively. 



III. BINDING ENERGY OF X c AND Xb 

The in-medium potential ([5]) resembles to the hydrogen 
atom problem. The solution of the Schrodinger equation 
gives the eigenvalues for the ground states and the first 
excited states in charmonium (J/ip, ip' etc.) and bot- 
tomonium (T, T' etc.) spectra : 



&n — K 

n z 



Er = 



rriQC 



(0) 



where mQ is the mass of the heavy quark and Ej is the 
energy of the QQ state in its first Bohr orbit. The allowed 
energy states for QQ are E n = — Ei, — j-, These ener- 
gies are known as the ionization potentials/binding ener- 
gies for the nth bound state. They become temperature- 
dependent through the Debye masses and decrease with 
the increase in temperature. 

Apart from the ground and the first excited states, 
there are other important states (lp) in the charmonium 
and bottomonium spectra viz Xc and Xb which contribute 
significantly in the suppression of ground state quarkonia 
(J/4> and T) in RHIC experiments through their decays 
into J/ip's and T's. Although both ip' and \ c are the 
first excited states of the charmonium spectra but they 
are not degenerate. In fact, ip' IS more massive than 
Xc but Xc is more tightly bound than ip'. So the entire 
binding energy of Xc will not come from the above cal- 
culation, the additional contribution will come from the 
spin-dependent quantum corrections. 

Some authors have studied the relativistic two-particle 
Coulomb problem, based on approximations to the 
Bether-Salpter equations. Others have started with ef- 
fective Lagrangians based on perturbativc expansions of 
the relativistic Lagrangians. However, we choose the 
variational methods [Hi, H3, H3| where coupled integral 
equations for a relativistic two-fermion system are de- 
rived variationally within the Hamiltonian formalism of 
quantum electrodynamics, using an improved ansatz that 
is sensitive to all terms in the Hamiltonian. The equa- 
tions are solved approximately to determine the eigen- 
values and eigenfunctions, at arbitrary coupling, for var- 
ious states of the two-particle system. In the variational 
treatment of the relativistic two-fermion bound-state sys- 
tem in QED, the total energy in a quantum state (nlsj) 
consists of Bohr like terms, relativistic correction in the 
kinetic energy and most importantly the spin-dependent 
terms which take into account the non-degeneracy be- 
tween the sub-states. The total energy up to fourth order 
in a is written for a hydrogen-like potential [32| 



Enlsj = 2mQ 



1 " 

-/i— + AK nl 



AV n 



(7) 
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where 



AK 



Til 



4 4 



mi \(2l + l)n 3 



3 
4r? 



(8) 



is the a 4 correction to the kinetic energy and the correc- 
tion to the spin-dependent potential energy is 



AV r 



nlsj 



9 nlsj ' 



(9) 



where the coefficients 6 nlaJ for the different quantum 
states {nlsj) are tabulated in Ref.[32||. We have taken 
their values for the Coulombic case only. The fine struc- 
ture constant (a) in QED will be replaced by the effective 
charge (2o/rn 2 D ) in our model. Using the appropriate 
values of the quantum numbers and the coefficients cor- 
responding to if)' and \c states in charmonium spectra 
(T' and Xb in bottomonium spectra), we obtain the cor- 
rection term which is to be added to the binding energy 
of ip' is 



pcorr 



m c ,ba 
96 



6m 8 D 



(10) 



So the binding energy of Xc(Xb) is 

E( Xc , Xb ) = EW,V)+El 



2 a 2 

!+o-4 
A m 



D 



— , (11) 



where mjj is the Debye mass for which we choose a gauge 
invariant, non-perturbative form [43|. Recently Kajantie 
et. al [43l ] obtained it by computing the non-perturbative 
contributions of 0{g 2 T) and 0{g 3 T) from a 3-D effective 
field theory as 



NP 



LO 1 N c g 2 T 
n c 9 T + d N N 



In 



L D 



g 2 T 
g 3 T + 0(g 4 T) 



(12) 



where the leading order (LO) perturbative result, 



g(T)T^j ^j- + has been known for a long 

time [iij]. The logarithmic part of the 0(g 2 ) correction 
can be extracted perturbatively [4{|, but cn c and the 
higher order corrections are non-perturbative. We wish 
to explore the effects of the different terms in the Debye 
mass on the binding energy of \c and Xb- We have used 
the two-loop expression for the QCD coupling constant at 
finite temperature from Ref. [4(| and the renormalization 
scale from Ref.[47j. 

The effects of each terms in the Debye mass (fT2"j) cannot 
always be explored separately due to the following rea- 
son: In the weak coupling regime, the soft scale (~ gT) 
at the leading-order related to the screening of electro- 
static fields is well separated from the ultra-soft scale 
(~ g T) related to the screening of magnctostatic fields. 
In such regime, it appears meaningful to sec the contri- 
bution of each terms in the Debye mass separately. But 



when the coupling becomes large enough (which is indeed 
the case), the two scales are no longer well separated. So 
while looking for the next-to-leading corrections to the 
leading-order result from the ultra-soft scale, it is not a 
wise idea to stop at the logarithmic term, since it be- 
comes crucial the number multiplying the factor 1/g to 
establish the correction to the LO result. In fact the De- 
bye mass in the NLO term is always smaller than the LO 
term because of the negative (logarithmic) contribution 
(log(l/<7)) to the leading-order term, while the full correc- 
tion (all g 2 T terms) to the Debye mass results positive. 
So, wc consider only three forms of the Debye masses, 
viz. leading-order result (m^-*), non-perturbative form 
(to^, p ), and lattice parametrized form (mj^ = 1.4m^°) 
to study the dissociation phenomena. 



Thus we have finally computed the binding energies 
for Xc and Xb and plotted them in Figs. 1 and 2, re- 
spectively where different curves denote the choice of the 
Debye masses used to calculate the binding energy from 
Eq. pTj) . We consider three cases for our analysis: pure 
gluonic, 2-flavor and 3-flavor QCD. There is a common 
observation in all figures that the binding energies show 
strong decrease with increase in temperature. In partic- 
ular, binding energies obtained from and mj^ give 
realistic variation with the temperature. The tempera- 
ture dependence of the binding energies show a quanti- 
tative agreement with the results based on the spectral 
function technique calculated in a potential model for 
the non-relativistic Green's function On the other 
hand, when we employ non-perturbative form of the De- 
bye mass (m NP ) the binding energies become unrealisti- 
cally small compared to its zero temperature value and 
also compared to the binding energies employing rrJff 



and 



D ' 



This anomaly can be understood by the fact 



that the value of rn^, p is significantly larger than both 
and mj^ so that the binding energies become sub- 
stantially smaller. This observation indicates that the 
present form of the non-perturbative corrections to the 
Debye mass may not be the complete one, the situation 
may change once the 0(g 4 T) non-perturbative contribu- 
tions to Debye mass are incorporated and then evaluate 
the binding energy. Thus, the study of temperature de- 
pendence of binding energy is poised to provide a wealth 
of information about the nature of dissociation of quarko- 
nium states in a thermal medium which will be reflected 
in their dissociation temperatures discussed in the next 
section. 



In addition, we take advantage of all the available lat- 
tice data, obtained not only in quenched QCD (TV/ = 0), 
but also including two and, more recently, three light fla- 
vors. We are then in a position to study also the flavor 
dependence of the dissociation process, a perspective not 
yet achieved by the parallel studies of the spectral func- 
tions, which are only available in quenched QCD. 
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FIG. 1: Dependence of Xc binding energy (in GeV) on temperature T/T c . 



IV. DISSOCIATION TEMPERATURES 

It has been customary to consider a state dissociated 
when its binding energy becomes zero. In principle, a 
state is dissociated when no peak structure is seen, but 
the widths shown in spectral functions from current po- 
tential model calculations are not physical. Broadening 
of states as the temperature increases is not included in 
any of these models. In [48[ authors have argued that no 
need to reach zero binding energy {E^- m = 0) to disso- 
ciate, but when £bin < T a state is weakly bound and 
thermal fluctuations can destroy it. However, others have 
set a more conservative condition for dissociation [4^ |: 
2E bin (T) < r(T), where T(T) is the thermal width of 
state. However, we now calculate the upper bound of 
the dissociation temperature by the condition for disso- 
ciation: \E Xc J wT, 

where the string tension (cr) is 0.184 GeV 2 and criti- 
cal temperatures (T c ) are taken as 270MeV, 203MeV 
and 197 MeV for pure gluonic, 2-flavor and 3-flavor QCD 
medium, respectively[50j. However, the choice of the 
mean thermal energy (T) is not rigid because even at low 
temperatures T < Td (say) the Bosc/Fermi distributions 
of partons will have a high energy tail with partons of me- 
chanical energy > \E Xc b |. The dissociation temperatures 
for \ c and Xb are listed in Table I with the Debye mass 
in the leading-order. It is found that Xc's are dissociated 
at 1.1T C , 1.31T C , and 1.26T C for the pure, 2-flavor, and 
3-flavor QCD, respectively whereas Xf>' s are dissociated 
comparatively much higher temperature which seems jus- 
tifiable. This is perhaps the first observation in the lit- 
erature on the flavor (system) dependence of the disso- 
ciation temperature. This dependence is essential while 



TABLE I: Upper(lower) bound on the dissociation 
temperature(TD) for Xc and Xb (in unit of T c ) using 
the leading-order term in Debye mass, m'ff '. 



State 


Pure QCD 


N f =2 


N f =3 




1.10 (0.89) 


1.31 (1.06) 


1.26 (1.02) 


Xb 


1.38 (1.10) 


1.64 (1.31) 


1.57 (1.26) 



TABLE II: Same as Table I but with the lattice parametrized 
form of the Debye mass m^. 



State 


Pure QCD 


N f = 2 


JV> = 3 


Xc 


0.77 (0.59) 


0.93 (0.73) 


0.90 (0.71) 


Xb 


0.98 (0.77) 


1.17 (0.93) 


1.14 (0.90) 



calculating the screening energy density (energy density 
at the dissociation temperature) in various descriptions 
of QGP (Nf = 0, 2, 3) for the study of J/tp survival in an 
expanding QGP. On the other hand, employing lattice 
parametrized form, mj^ we obtain the values (Table II) 
much smaller than the leading-order results where Xc is 
dissociated below T c and Xb is dissociated just above T c . 
At last, when we use non-perturbative form of the Debye 
mass, the dissociation temperatures come out to be un- 
realistically small. Summarizing the results, we conclude 
that as we move from perturbative to non-perturbative 
domain, the binding energies are becoming smaller and 
smaller. As a result the dissociation temperatures ob- 
tained are also becoming smaller. This is due to the 
hierarchy in the Debye masses: m^f* < <C rrJj^ . In 

fact, mJfy is 1.4 times greater than while rn^ is 

much greater than both mjy and mj^. 

However, if we treat the partons in high temperature to 
be relativistic, we could replace the mean thermal energy 
by 3T (instead of T) to obtain the lower bound for the 
dissociation temperatures. It is found that all entries in 



G 




FIG. 2: Dependence of Xf> binding energy (in GeV) on T/T c . 



TABLE III: Upper bound of the dissociation temperatures 
with T c = 192 MeV and the lattice parametrized form of the 
Debye mass [48l ]. 



State 


Xc 


Xb 




< 0.9T C 


1.2T C 



Table I and II have been decreased by 30% approximately 
giving the lower bound of the dissociation temperatures 
(inside the first bracket). To compare our results quanti- 
tatively with the recent results [48[ based on the spectral 
function technique calculated in a potential model with a 
similar description of the system (for 3-flavor QCD with 
T c =192 MeV), we tabulated the upper limit on the disso- 
ciation tem per atures with the same form of Debye mass 
used in Ref . [481 ] in Table III giving a good agreement with 
their results. 

Finally, it is learnt that inclusion of non-perturbativc 
corrections to the Debye mass (ni^ p ) leads to unusually 
smaller value of the dissociation temperatures for both 
Xc and Xb- This does not immediately imply that the 
non-perturbative effects should be ignored. It is rather 
interesting to investigate the disagreement between the 
non perturbative result obtained with a dimensional- 
reduction strategy and the Debye mass arising from the 
Polyakov-loop correlators. Only future investigation may 
throw more light on this issue. 

V. CONCLUSIONS 

In conclusion, we have studied the dissociation of lp 
states in the charmonium and bottomonium spcctra(x c 
and Xfc) m the hot QCD medium. We have employed 
the medium modified form of the heavy quark-potential 
in which the medium modification causes the dynam- 



ical screening of color charge which, in turn, leads to 
the temperature dependent binding energy of ip' and T\ 
We have then studied the temperature dependence of the 
binding energy of the x c and X b states in the pure gauge 
and realistic QCD medium by incorporating the fourth- 
order corrections (in the screened charge a e ^ = 2<j/m 2 D ) 
coming from the spin dependent terms to the binding 
energies of tp' and T' states, respectively. For this pur- 
pose, we have adopted a formulation [13] , in which a vari- 
ational treatment of the relativistic two-fermion bound- 
state system in QED (32j has been developed to compute 
the spin-dependent corrections. 

Next we have determined the dissociation tempera- 
tures employing the Debye mass in leading-order and the 
lattice parametrized form. Our estimates are consistent 
with the findin g of recent theoretical works based on po- 
tential models [48|]. We have further shown that inclu- 
sion of non-perturbative contributions to the Debye mass 
lower the dissociation temperatures substantially which 
looks unfeasible to compare to the spectral analysis of lat- 
tice temporal correlator of mesonic current. This leaves 
an open problem of the agreement between these two 
kind of approaches. This could be partially due to the 
arbitrariness in the criteria/definition of the dissociation 
temperature. To examine this point we have estimated 
both the upper and lower bound on the dissociation tem- 
peratures by fixing the mean thermal energy T and 3T, 
respectively. Thus, this study provides us a handle to 
decipher the extent upto which non-perturbative effects 
should be incorporated into the Debye mass. 

In brief, we obtained the analytic forms for the binding 
energies and the dissociation temperatures of Xc and Xb- 
This enable us to investigate their flavor dependence and 
temperature dependence. We have estimated the upper 
bound on the dissociation temperatures of Xc and Xb- We 
found that these estimates obtained by employing the 
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lattice parametrized Debye mass show good agreement 
with the prediction in [48|. On the other hand, these 
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